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ABSTRACT 



We investigate the angular correlation function, oj(9), and third and fourth-order 
hierarchical moments of a large sample of ~ 70000 galaxies of apparent magnitude 
18.5 < R < 23.5. The data consists of 47 red-band INT Wide Field Camera CCD 
\ images, forming two widely separated fields of total areas 1.01 deg 2 and 0.74 deg 2 . 

Galaxy clustering is detected with a high significance of ~ 10c Over the R = 21 to 
R = 23.5 range of magnitude limits, the angular correlation functions approximately 
follow 6>~ - 8 power-laws at 9 > 5 arcsec, with amplitudes consistent with previous 
surveys and best-fitted by a luminosity evolution model in which galaxy clustering is 
approximately stable in proper co-ordinates (e = 0). Assuming the redshift distribution 
from our pure luminosity evolution model and the present-day galaxy correlation radii 
• from Loveday et al. (1995), we estimate the clustering evolution as e = 0.021q; 3 * for 

OO . bth fields combined or e = — 0.29^q'27 for the larger of the two fields which is thought 

ON ' to be the better in data quality. 

^ (— I On the larger of our fields, lu(9) at 2 < 9 < 5 significantly exceeds a uj(9) cx 9~ o s 

O power-law of the amplitude measured at 9 > 5 arcsec. If this excess of close pairs is 

I ' due to interacting and merging galaxies, we estimate that it is consistent with the local 

fraction of galaxies in close (< 19/i _1 kpc) pairs combined with a merger/interaction 
^ ' rate evolving as (1 + z) m with m = 2.01:^52 

\ We derive the hierarchical moments s-sQi) and s^iG) from the counts- in-cells of the 

18.5 < R < 23.5 galaxies. We find relatively steep slopes of approximately 9~ 0A for 
S3 (6*) and 9~ 6 for S4(9), similar to those of the hierarchical moments of 17 < B < 20 
galaxies in the EDSGCS (Szapudi et al. 1996) at comparable physical scales. The 
■ s 3($) an d £4(0) of the 18.5 < R < 23.5 galaxies show only a small reduction in 

normalization relative to those of the less deep EDSGCS survey, consistent with the 
expectations from A^-body simulations of the evolving mass distribution. This indicates 
there is little change in the mean linear or non-linear biasing of galaxies from z ~ 
to z ~ 1, and hence appears to favour luminosity evolution for galaxies over 'transient 
starburst dwarf models, and supports the interpretation of the stable galaxy clustering 
as evidence for a low ft Universe. 

Key words: 

surveys - galaxies: clusters: general - galaxies: evolution - large-scale structure of 
Universe 



1 INTRODUCTION 

The clustering of galaxies on the sky is most often described 
in terms of the angular correlation function, w(6 l ), a projec- 
tion onto the sky plane of the three-dimensional two-point 
correlation function £(r). For galaxies of separation r in 
proper co-ordinates, at redshift z, the two-point correlation 
function is often parameterized as 



ar,z) = (^y(l + z)-^ (1) 



where ro is the correlation radius, 7 the slope, and e pa- 
rameterizes the evolution of clustering (clustering stable in 
proper co-ordinates being e = 0). The angular correlation 
function will be a power-law of slope g - ^ -1 ' and amplitude 
given by the Limber's formula (e.g. Phillipps et al. 1978) 
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integration of £(r, z) over the galaxy redshift distribution 
N(z). 

Photographic galaxy surveys indicated that u>{9) ap- 
proximately follows a #~ 0,8 power-law, with a break at very 
large scales, and an amplitude which decreases with increas- 
ing survey depth due to the effects of projection (e.g. Groth 
and Peebles 1977; Maddox et al. 1990). During the early 
1990s, a number of CCD surveys indicated that the u>(6) am- 
plitude of faint (B > 23.5) galaxies falls well below the scal- 
ing expected for stable clustering and a galaxy N(z) given 
by a non-evolving model (e.g. Efstathiou et al. 1991; Roche 
ct al. 1993). This was interpreted (Roche et al. 1993, 1996) 
as evidence that N(z) at these magnitudes is more extended 
than the non-evolving prediction, implying an evolutionary 
brightening of L ~ L* galaxies at higher redshifts. Alterna- 
tively, it could mean that the clustering of galaxies evolves 
rapidly (i.e. e > 1) and so was much weaker at even mod- 
erate redshifts. Another possibility was that that the faint 
blue galaxies at B > 23 were a new population of blue dwarf 
galaxies at moderate redshifts, with much weaker clustering 
than L ~ L* spirals (Babul and Rees 1992; Efstathiou 1995). 

It was then found that the redder galaxies at these mag- 
nitudes gave a much (~ 0.5 dex) higher to(6) amplitude 
than the bluer galaxies (Neuschaefer et al. 1995; Roche et al. 
1996). This would argue against a strong evolution of galaxy 
clustering as an explanation of the low u>{ff) amplitudes, but 
would be consistent with L* evolution models, where the ex- 
tended tail of N(z) would consist of blue star-forming giant 
galaxies. However, it would also be consistent with intrinsi- 
cally weakly clustered blue dwarf galaxies at lower redshift. 

Soon afterwards, the importance of L* evolution was 
established directly when deeper redshift surveys with the 
Keck telescope (e.g. Cowie, Songaila and Hu 1996; Steidel et 
al. 1996) and Hubble Deep Field photometry (e.g. Madau et 
al. 1996; Bershady et al. 1997) identified many galaxies at 
z ~ 1-3 and confirmed there is a brightening of the galaxy 
L* at these redshifts. 

As the luminosity evolution is becoming better deter- 
mined by other data, the emphasis of studying faint galaxy 
clustering is changing from the use of u>{ff) as a probe of 
N(z) towards a more detailed investigation of the evolution 
of clustering and large-scale structure, the clustering prop- 
erties of different types of galaxy, and the possible effects 
of interactions between galaxies. In particular, there is still 
much to learn about the rate of clustering evolution, e, which 
is sensitive to the cosmological parameters and also to any 
biasing of the luminosity distribution relative to the under- 
lying mass distribution. Models (e.g. Colin, Carlberg and 
Couchman 1997) and direct observations from redshift sur- 
veys (e.g. Le Fevre et al. 1996; Shepherd, Carlberg and Yee 
1997; Carlberg et al. 1997) generally span a wide range of at 
least < e < 1.5, although models in which the bias factor 
increases rapidly with redshift (e.g. Matarrese et al. 1997) 
may predict galaxy clustering evolution closer to a comoving 
model with e = 7 — 3 ~ —1.2. 

The first deep CCD surveys covered small (< 0.1 deg 2 ) 
areas of sky, on which the detection of galaxy clustering 
was only ~ 3a, but with large format CCDs, such as those 
in the Wide Field Camera installed on the Isaac Newton 
Telescope in 1997, it becomes possible to survey > 1 deg 
areas and study faint galaxy clustering in more detail. In this 
paper we present a study of galaxy clustering to a limit R — 



23.5 on two fields with total areas 1.01 deg 2 and 0.74 deg 2 . 
The larger fields will greatly reduce the 'integral constraint' 
correction (see Section 4), and the increase in sample size 
will provide better constraints on the clustering evolution. 

Previous, less deep, surveys suggest that > 1 deg 2 areas 
can provide more than just the ui(9) amplitude. For example, 
Infante and Pritchet (1995), using a 2deg 2 mosaic of photo- 
graphic plates, reported some flattening of the slope of uj(6) 
between R ~ 21 and their limit of R ~ 23, while Infante, de 
Mello and Mentaneau (1996), using a similar area made up 
of short-exposure CCD images, found the R < 21.5 galaxy 
lo(8) to be significantly higher in amplitude at 2 < 6 < 6 
arcsec than at larger scales. With a large area, it may also be 
possible to measure the higher-order correlations functions, 
which are related to the higher-order moments (skewness, 
kurtosis etc.) of the galaxy counts-in-cells, in the same way 
that u>{0) is related to the variance. These statistics provide 
additional information on the nature of galaxy clustering but 
until now have been studied only at much shallower limits 
(Gaztanaga 1994; Szapudi, Meiksin and Nichol 1996). 

We assume Ho = 50 km s _1 Mpc _1 and go = 0.05 
throughout, but give some quantities in units of h = Ho/WO 
km s~ Mpc - . In Section 2 of this paper we describe the 
observational data and its reduction and calibration, and in 
Section 3 the number counts of the galaxies and stars. Sec- 
tion 4 describes the analysis of the galaxy lo(0), and Section 
5 presents the results. Section 6 describes models of galaxy 
evolution and clustering and compares these with our obser- 
vations and other data. In Section 7 we discuss lu(0) at small 
scales of a few arcsec and its use in investigating merger rate 
evolution. In Section 8 we investigate the higher-order mo- 
ments of the sample and compare with results from less deep 
surveys. Section 9 discusses in detail the implications of all 
these results. 



2 DATA 

2.1 Observations 

Our data were acquired on the nights of 4th and 5th June 
1997, on a service run of the recently installed Wide Field 
Camera (WFC), at the prime focus of the 2.5 metre Isaac 
Newton Telescope on La Palma. At that time the WFC was 
fitted with a 2 x 2 array of Loral CCDs, each with 2048 x 2048 
pixels. In this configuration the pixelsize corresponds to 0.37 
arcsec. The WFC should then have provided a total area of 
approximately 638 arcmin 2 with each pointing, but one of 
the four CCD chips was nonfunctional and another suffered 
a severe loss of sensitivity towards the edges, reducing the 
usable area to ~ 440 arcmin 2 . 

The WFC imaged 21 positions, 12 arranged in a 3 x 4 
grid referred to as field 'e', and 9 in a 3 x 3 grid referred 
to field 'P. Field 'e' is centred at R.A. 14:00 hours, declina- 
tion zero, which is galactic longitude 337.01° and latitude 
+58.26°, and field T is centred at R.A. 21:00 hours, decli- 
nation — 10°, which is galactic longitude 38.71° and latitude 
—33.06°. Each of the 21 pointings was exposed twice for 
900 seconds, with the WFC Harris Red filter. For photo- 
metric calibration, 5 second exposures were taken of several 
fields containing standard stars from the catalog of Landolt 
(1992). 
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Table 1. Positions (equinox 2000.0) of the 16 pointings of the 
WFC and estimates of the resolution on the images 



Field name 


H.A. 


Dec. 


r WnM (arcsec ) 


el 


"i a nn nn nn 

14:00:00.00 


i / 1 ( i / 1 j i 1,1, • i 
+00:00:00.3 


1.43 


e2 


13:59:59.98 


+00:25:39.9 


1.41 


e4 


14:00:00.00 


-00:51:19.7 


1.51 


e7 


14:01:38.05 


+00:00:00.3 


1.70 


e8 


14:01:37.94 


-00:24:29.7 


1.73 


e9 


14:01:37.84 


-00:49:01.3 


1.62 


elO 


1 3-58-22 1 


-00-00-00 1 


1.51 


ell 


13:58:22.07 


+00:24:29.3 


1.63 


el2 


14:01:38.16 


+00:24:29.3 


1.88 


fl 


21:00:00.09 


-10:00:00.2 


1.92 


f2 


20:59:59.97 


-10:24:29.9 


1.89 


f4 


21:01:39.00 


-09:35:29.8 


1.72 


f6 


21:01:38.98 


-10:24:25.7 


2.19 


f7 


20:58:21.07 


-09:35:29.2 


1.68 


f8 


20:58:21.01 


-09:59:58.5 


1.55 


f9 


20:58:21.03 


-10:24:28.7 


1.78 



Some further sections of the data were rendered unus- 
able by autoguider malfunctions and a period of poor seeing. 
We were able to use a total of 47 CCD frames (27 in field 
'e' and 20 in field 'f ') from 16 pointings of the WFC at co- 
ordinates given in Table 1. 



2.2 Reduction and source detection 

The initial processing of the data was carried out using IRAF. 
For each of the three functional CCD chips, a bias frame 
was produced and subtracted from all images. One of the 
three CCD chips showed a significant falloff in sensitivity 
towards the edges. We trimmed these images approximately 
where the sensitivity fell to 0.64 of its central value, reducing 
the area by ~ 25%. For each chip, median flat fields were 
generated from the data using IRAF 'ccdcombine', by taking 
the median of all 32 exposures on the positions in Table 1, 
rejecting those pixels containing bright objects. All images 
were then divided by the respective normalized flat fields. 

The two exposures at each position were added using 
IRAF 'ccdcombine', using the 'crrej' option (rejecting the 
higher of the two values if discrepant by > 3a) to remove the 
majority of cosmic rays. Sky subtraction was performed by 
fitting high-order spline surfaces to the sky background on 
each image, using KAPPA 'surfit', and subtracting these from 
the data. The mean sky background intensity (using the 
calibration in Section 2.3) was 20.77 ±0.02 R mag arcsec -2 
on field V and 20.70 ± 0.02 R mag arcsec" 2 on field 'P. 

The Starlink PISA (Position, Intensity and Shape Anal- 
ysis) package, developed by M. Irwin, P. Draper and N. 
Eaton, was used to detect and catalog the objects, with 
the chosen detection criterion that a source must exceed an 
intensity threshold of 1.5a above the background noise (a 
being separately determined by PISA for each image) in at 
least 6 connected pixels. The mean detection threshold was 
25.72 ± 0.01 R mag arcsec -2 on field 'e' and 25.63 ± 0.02 R 
mag arcsec -2 on field 'P. PISA was run with deblend (ob- 
jects connected at the detection threshold may be split if 
they are separate at a higher threshold) and the 'total mag- 
nitudes' option, which estimates total intensity above the 




IB 20 
R magnitude 

Figure 1. The total/peak intensity parameter used for star- 
galaxy separation, plotted against R magnitude, for detections 
on one CCD frame within field 'e', showing their classification as 
stars or galaxies. The rise in the stellar locus at R > 18 is the 
result of saturation of the CCD chips imposing a upper limit on 
the measured peak intensity. 



background by fitting an elliptical aperture to each individ- 
ual detection and performing a 'curve-of-growth' analysis. 

We excluded data from circular 'holes' (radius 15 to 
100 arcsec) around many of the bright, saturated stellar im- 
ages in order to remove spurious noise detections. Radial 
profiles were fitted to several non-saturated stars on each 
image using 'pisafit', giving the FWHM estimates in Table 
1 . Evidently, the seeing was less than ideal during the period 
of observation, resulting in a FWHM averaging 1.60 ± 0.04 
arcsec on field 'e' and 1.82 ± 0.07 arcsec on field 'P, with 
much variation between the individual images. 

Star-galaxy separation was performed using plots of 
central against total intensity, normalized to the ratio from 
the fitted stellar profile. Figures 1 and 2 show the star-galaxy 
separation plots for one CCD frame from each of the two 
fields. The stellar locus remains separate from the galaxies 
only to R ~ 21-21.5; so detections fainter that this are all 
classified as galaxies. Even at 20.5 < R < 21.5 the reliability 
of star-galaxy separation appears to vary between images. 
Hence we only assume the separation is reliable to R = 20.5, 
and later (Sections 3 and 5) apply corrections to our results 
for the effects of star contamination at R > 20.5. 



2.3 Calibration 

For photometric calibration, the 900s exposures were inter- 
spersed with short exposures of three fields containing stan- 
dard stars (about 5 per CCD chip per field) from the catalog 
of Landolt (1992). After flat-fielding the standard star im- 
ages, the intensities of the standard stars were measured 
using IRAF aperture photometry. 

For each observation of a standard star, the photon 
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Figure 2. As Figure 1 for one image within field 'P. 



count rate C will relate to the catalogued red magnitude 
R and colour V — R as 



2.5 log(C) = R + n + r 2 X + r 3 (V - R) 



(2) 



where n is the zero-point outside the atmosphere, r 2 the 
atmospheric extinction at the zenith, X the airmass of each 
observation is and rz a colour-correction term expressing the 
difference between our passband and the Cousins R band of 
Landolt (1992). As the standard star observations were too 
few in number and covered too small a range of airmass to 
accurately measure ri, we estimated r 2 = 0.062 by integrat- 
ing a tabulated measurement of the extinction curve at La 
Palma over the R passband. 

The three CCD chips were calibrated separately. The 
colour term r$ was consistent with zero within the statistical 
errors, so the standard star observations were fitted with 



2.5 log(C) = R + n + 0.062X 



(3) 



giving n = -25.185±0.012, -25.135±0.005 and -25.369± 
0.004 for the three CCD chips. The photon counts for each 
detected source could then be converted to R magnitudes at 
the zenith with 



R = -2.5 log(C) - (n + 0.062X) 



(4) 



with n and X appropriate to the chip and airmass of ob- 
servation. 

A further correction is applied for Galactic dust extinc- 
tion. The Caltech NED database gives B-band extinctions 
of 0.10 and 0.19 magnitudes for objects near the centres of 
fields 'e' and 'f respectively, approximately in proportion 
to the Galactic HI column densities. Using the extinction 
curve of Mathis (1990), we estimate Ar — 0.65As, and ap- 
ply corrections of A(R) = -0.065 and A(R) = -0.124 to 
magnitudes from fields 'e' and 'f respectively. 



2.4 Astrometry 

Astrometric transforms were derived using of the APM 
database at the Institute of Astronomy in Cambridge. Field 
'e' lay within the measured UKST plate F865 and field 'f 
within the measured UKST plate F742. We identified 12 cat- 
alogued stars on each of the 47 CCD images, and using their 
R.A. and Dec. fitted astrometric transforms (with residuals 
typically only ~ 0.5 arcsec) using the I RAF 'pltsol' routine. 
The pixel co-ordinates of all detections were then converted 
to R.A. and Dec. The catalogs of detections on the 47 indi- 
vidual images could then be combined into two composite 
catalogs for the 'e' and 'f fields. 

Most of the 16 exposures overlap with adjacent point- 
ings by a few arcminutes. For the areas of overlap, we 
adopted the positions and magnitudes derived from the im- 
age with the better seeing (Table 1) when combining the 
detection catalogs. Figures 3 and 4 show the R.A and Dec 
positions of the brighter galaxies in the composite catalogs 
for fields 'e' and 'f. In calculating angular separations be- 
tween galaxies we assume A(#x)/A(R.A) to be unity on 
field 'e' and cos(10°) on field 'f, neglecting the very small 
(< 0.02%) effect of sky curvature within each field. After 
taking into account the overlap regions and holes, the areas 
covered by usable data amount to 1.012 deg 2 for field 'e' and 
0.740 deg 2 for field 'P. 



2.5 Photometric Matching in Overlap Areas 

The CCD frame overlap areas were made use of in two ways 
- firstly, to estimate the photometric errors for individual 
faint galaxies by matching the two detections of the same 
objects and comparing the magnitudes. The mean scatter 
between the two detections gave errors as a{R) = 0.14 mag 
for R < 22 detections, 0.22 mag at 22 < R < 23 and 0.24 
mag at 23 < R < 24. 

Secondly, by calculating a mean offset between the mag- 
nitudes of the same objects in each overlap, we can check 
for photometric shifts between the exposured, typically to 
~ 0.05 mag accuracy. Significant offsets were seen for three 
exposures only - e7, e8 and f6 were estimated to have suf- 
fered losses of 0.18±0.06 mag, 0.41±0.09 mag and 0.09±0.04 
mag respectively, in comparison to the rest of the data, pre- 
sumably as a result of poor atmospheric conditions. Ap- 
propriate corrections (—0.18, —0.41 and —0.09 mag) were 
therefore applied to all magnitudes from these fields. The 
overlap areas will reveal only the worst of the photometric 
offsets, and, inevitably, there will be some remaining scatter 
in the exposure zero-points - this is investigated further in 
Section 5.3. 



3 GALAXY AND STAR COUNTS 

Figure 5 shows the differential number counts of objects 
classed as galaxies on the 'e' and 'P fields (error bars are 
derived from the scatter in sufface density between the in- 
dividual CCD frames making up each field), together with 
galaxy counts from some previous CCD surveys in similar 
passbands, and the PLE and non-evolving models described 
in Section 6. Comparison with deeper galaxy counts suggests 
our detection of galaxies is virtually complete to R = 23.5, 
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Figure 5. Differential galaxy number counts, in 0.5 mag intervals 
of R magnitude, for our two survey fields compared to previous 
CCD surveys in very similar passbands and the predictions of our 
PLE (dashed) and non-evolving (dotted) models. 



Figure 3. The positions of detections classed as 18.5 < R < 22.5 
galaxies on the 27 CCD frames comprising field 'e' 
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Figure 4. The positions of detections classed as 18.5 < R < 22.5 
galaxies on the 20 CCD frames comprising field 'f 



but becomes significantly (~ 30% on field 'e', ~ 45% on 
field 'f') incomplete at 23.5 < R < 24, before our counts 
turn over at R > 24. As incompleteness can unpredictably 
affect measures of clustering, our analysis is limited to to 
the R < 23.5 galaxies. 

At 18.5 < R < 23.5, the field 'e' galaxy counts are close 
to both previous observations and the PLE model, whereas 
the field 'f counts are significantly higher at 20.5 < R < 
22.5. The galaxy counts also show an excess above the mod- 
els at bright magnitudes of R > 18.5, on both fields but 
much more so on field 'P. These discrepancies are the result 
of star-contamination. 

We estimate the contamination by stars faintward of 
the limits of star-galaxy separation (R ~ 20.5-21) using the 
Galactic star count models from Reid et al. (1996). The R- 
band models for Galactic latitudes b = 60° and b — 45° are 
adopted for fields 'e' and 'f ' respectively and normalized to 
fit our observed star counts (figure 6) at 18.0 < R < 20.5. 
The fitted normalizations corresponded to 1087 ± 76 stars 
mag _1 deg -2 on field 'e' and 2791 ± 136 stars mag _1 deg -2 
on the lower Galactic latitude field 'P, with errors from the 
scatter in star counts between the CCD frames within each 
field. 

We then (i) assume stars are classified correctly to 
R = 20.5, (ii) estimate the number of R > 20.5 stars by 
summing the normalized model star-counts from R = 20.5 
to the magnitude limit being considered, and (iii) subtract 
the number of R > 20.5 detections already classified as stars. 
To R — 23.5, this gives the total number of contaminating 
stars as 4306 ± 301 on field 'e' and 8433 ± 411 on field 'P. 

Brightward of R ~ 18-18.5 the star counts from both 
fields fall below the plotted models, approximately mirror- 
ing the excess seen in the galaxy counts. Saturation of the 
brighter star images in our data evidently results in sig- 
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Figure 6. Differential star number counts in 0.5 mag intervals 
of R magnitude, with star count models (dotted line for field 'e', 
dashed line for field 'f') from Reid et al. (1996), normalized to the 
observed count at 18 < R < 20.5. 



nificant numbers becoming misclassificd as galaxies. As we 
are interested in the clustering of much fainter galaxies, we 
remove this problem simply by excluding all R > 18.5 detec- 
tions from our analysis. This leaves a total of 80293 objects 
in our dataset classed as 18.5 < R < 23.5 galaxies, of which 
we estimate 12739 ± 509 are stars, leaving 67554 ± 509 gen- 
uine faint galaxies. 



4 CALCULATING w(0) 

The angular correlation function u(6) was calculated sepa- 
rately for fields 'e' and 'f. For field 'e' we used all galaxies 
brighter than a series of limits from R = 21.0 to R = 23.5, 
and on field 'f used the same faint limits but also excluded 
R < 19.0 objects. To calculate uj{9) for a field containing 
Ng a i galaxies in the chosen magnitude range, the separa- 
tions of all \Ng(N g — 1) possible galaxy-galaxy pairs were 
counted in bins of width A(log 9) = 0.2, giving the function 

For each field N r — 100000 points were placed at ran- 
dom over the area covered by real data. i.e. avoiding any 
gaps or holes. The separations of the N r N g galaxy-random 
pairs, taking the real galaxies as the centres, were similarly 
binned, giving N rg (6i) and likewise the separations of the 
^N r (N r — 1) random-random pairs were counted to give 
N rr {9i). 

Defining DD — N gg (9i), and DR and RR as the galaxy- 
random and random-random counts normalized to have the 



same summation over all 9 as DR, i.e. DR = 



N, 



gr{9 l ) 



and RR 



Ng(Ng-l) 



N rr (9i), we might then estimate co(9) 



N R (N r -l - 

for each bin simply as (Roche et al. 1993, 1996) 



"(*0 = TTd 



DD 
~DR 



(5) 



However, Landy and Szalay (1993) presented a new estima 
tor of to(9), said to give the smallest possible statistical er- 
rors, and to optimally remove edge effects (Bernstein 1994) 
of the form 

DD — 2DR + RR 



RR 



(6) 



Ratcliffe et al. (1998) compared these clustering estimators 
using simulated galaxy catalogs, and confirmed that the 
Landy and Szalay (1993) estimator was significantly more 
accurate than that in equation (5) for the correlation func- 
tions at larger separations, so it will be used for the analysis 
of this paper. 

Figures 7 and 8 show u)(6) calculated with equation (6) 
for fields 'e' and 'f respectively. The error bars in u>(9) esti- 
mates can be much larger than the Poisson error (DD -0,5 ) 
and depend in a complex fashion on the field geometry, 
galaxy surface density and higher-order correlation func- 
tions. We estimate the error bars using a 'jackknife' method, 
in which u>(9) is recalculated 10 times for each field, each 
time excluding just one of the CCD frames (randomly cho- 
sen) from the analysis. As these subsamples will be smaller 
than the full fields by only one part in 27 (field 'e') or 20 
(field 'f'), the effects of altering the field geometry will be 
minimized. For field 'e', the error on u(9) for the full field of 
27 images was estimated by multiplying the scatter in u>{9) 
between these 10 subsamples by \/26 (as only one part in 
26 of the data in changed from one subsample to the next) 
x \f^j (as the full field contains 27 rather than 26 images) 
= 5.00. For field 'f, the error for the full field of 20 images 
was similarly estimated by multiplying the scatter between 
the subsample u>(9) images by \/l9 x \fl^ = 4.25. 

If the real galaxy lj(6) is of the form A9~ ' 8 , the ob- 
served tu(9) will follow the form = A(9~ 08 - C), with 
amplitude A (defined here at a one-degree separation), and 
a negative offset AC known as the integral constraint. This 
offset results from the restricted area of the observation, and 
can be estimated by doubly integrating an assumed true 
u)(6) over the field area fi, i.e. 



AC= T? 



uj{9)dn 1 dQ 2 



(7) 



This calculation can be done numerically using the random- 
random correlation. Assuming u>(9) — A9~ ihS , 



C 



^A rr (fl)fl- -' 

ZN rr {9) 



We then fit u(9) with a function w(0) = A(9~ 



(8) 
-C), 

but find that, at least on field 'e', the two lo(9) points at 
2 < 9 < 5 arcsec lie significantly above a 9~ S power-law 
fitted at larger scales. Infante et al. (1996) had previously 
reported that the u(8) of a large sample of R < 21.5 galaxies 
appeared to retain a 9~ ' 8 slope at 2 < 8 < 6 arcsec but was 
increased in amplitude by a factor 1.8 relative to the u>(9) 
at 9 > 6 arcsec. This result, indicating an excess of close 
pairs, was interpreted as the result of galaxy interactions 
and mergers. 

In order to quantify the effects of interactions in our 
deeper survey, and at the same time minimize their influence 
on the measurement of larger scale galaxy clustering, we fit 
our u>(9) with two 9~ ' 8 power-laws, one at 69 > 5 arcsec to 
derive an amplitude A, another at 2 < 9 < 5 arcsec to derive 
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Table 2. Observed ui(6) amplitudes at 9 > 5 arcscc (A) and at 
2 < 9 < 5 arcsec (A s ), in units of 10~ 4 at 1°, for the galaxies, 
numbering N ga i, with R magnitudes from R = 18.5 to a series of 
faint limits. 



R limit 


A 


A a 


Ngal 






Field 'e' 




21.0 


38.93 ± 4.97 


55.44 ±5.57 


4395 


21.5 


27.00 ±2.00 


40.15 ±5.58 


7643 


22.0 


20.20 ± 1.32 


38.61 ±3.31 


12665 


22.5 


17.66 ± 1.54 


33.31 ±2.94 


19890 


23.0 


15.73 ± 1.85 


27.28 ± 2.49 


30284 


23.5 


16.07 ± 2.08 


23.74 ± 1.89 


44307 






Field 'f 




21.0 


26.52 ±2.77 


27.50 ±6.71 


5274 


21.5 


15.29 ±2.59 


12.97 ±2.77 


8877 


22.0 


13.06 ± 1.81 


11.97 ±2.36 


13567 


22.5 


12.00 ±2.11 


13.70 ± 1.87 


19415 


23.0 


12.59 ±2.18 


14.86 ±2.41 


26958 


23.5 


13.83 ±2.54 


16.60 ±2.50 


35986 



Table 3. Estimated fraction of stars f s contaminating the galaxy 
samples at each magnitude limit, and the u>(9) amplitudes at 
9 > 5 arcsec (A) and 2 < 9 < 5 arcsec (A s ) multiplied by a 
(1 — f s )~ 2 correction for star contamination. 



R limit 


fs 


A (corrected) 


A s (corrected) 






r lela c 




21.0 


U.035 / ± 0.U0Z5 


A 1 07 _l_ c: Q c; 

41.8 / ± o.oo 


£n an A- a c\c\ 
59. OZ ± 6.0U 


21.5 


0.1023 ± 0.0715 


33.50 ± 2.54 


49.82 ± 6.97 


22.0 


0.1217 ±0.0851 


26.19 ± 1.78 


50.05 ±4.40 


22.5 


0.1203 ± 0.0841 


22.82 ± 2.04 


43.04 ±3.89 


23.0 


0.1099 ± 0.0768 


19.85 ± 2.36 


34.43 ±3.20 


23.5 


0.0972 ± 0.0679 


19.72 ± 2.57 


29.13 ±2.36 






Field T 




21.0 


0.1342 ± 0.0066 


35.38 ±3.74 


36.69 ±8.97 


21.5 


0.2132 ± 0.0104 


24.70 ±4.23 


20.95 ±4.51 


22.0 


0.2446 ± 0.0120 


22.89 ±3.26 


20.98 ±4.20 


22.5 


0.2515 ± 0.0123 


21.42 ±3.83 


24.45 ±3.44 


23.0 


0.2447 ± 0.0120 


22.07 ±3.89 


26.05 ±4.30 


23.5 


0.2343 ± 0.0115 


23.59 ± 4.39 


28.31 ±4.35 



a 'small scale amplitude' A s (at 9 < 2 arcsec many pairs will 
be unresolved so the u>(9) cannot be used). 

The integral constraint then becomes ACi+A a C2 where 
Ci and C2 are, respectively, the summations of equation (9) 
at all angles 9 > 5 arcsec and at 2 < 9 < 5 arcsec only. 
For field 'e', Ci = 1.799 and C 2 = 0.00441 and for field 'f, 
Ci = 2.024 and C2 = 0.00517. As the integral constraint is 
greatly dominated by the large-scale u>(9), so we can neglect 
the contribution from any small-scale excess and assume the 
integral constraint to be A(Ci + C2). We estimate A by least- 
squares fitting A(9~ 0,li — (C1+C2)) to the observed ui(9) from 
5 arcsec to 7962 arcsec on field 'e' and 5024 arcsec on field 'f '. 
Once A has been determined, giving the integral constraint, 
A s is determined by fitting A s 9~ - 8 - A(Ci +C 2 ), to the two 
bins at 2 < 9 < 5 arcsec. 

The fitted functions, plotted on Figures 7 and 8, fit 
the field 'e' uj(9) at all separations, but the field 'f io{9) 
shows some evidence of an excess above the power-law at 
log 9 > —1.2 - we discuss a possible explanation in Section 
5.3. 

For comparison, we also calculate 'individual frame' 
u>{9), in which each CCD frame is treated as an indepen- 
dent field for which to(9) calculated using equation n above 
with N r — 20000, and the results averaged for the 27 frames 
of field 'e' and the 20 frames of field 'f. The integral con- 
straints are larger for the small fields of invididual frames, 
averaging C = 8.28. In Section 5.3 these estimates are com- 
pared with the full- field u>{9). 



timates due to the effects of star contamination. A fraction 
of (randomly distributed) stars f s contaminating the galaxy 
sample will dilute the observed lo(9) at all angles by a fac- 
tor of (1 — /s) -2 . Table 3 gives estimates of f s in our data 
from the fitted star count models of Section 3.4, and co(9) 
amplitudes A and A s corrected for star contamination by 
multiplying by (1 — fs)~ 2 , with errors which include the 
(small) contribution from the uncertainty in f s , added in 
quadrature to the ui(9) errors. 

These results show that, firstly, our detection of cluster- 
ing is of high significance, e.g. at R < 23.0 an estimated 8.5ct 
on field 'e' and 5.8<r on field 'f, or 10.3cr overall. Secondly, 
the two fields give reasonably consistent u>(9) amplitudes at 
9 > 5 arcsec, when the different star contamination is taken 
into account. Thirdly, at 2 < 9 < 5 arcsec the u>(9) am- 
plitude from fields 'e' is significantly higher than at larger 
separations, by > 3.5a at R = 22-23 limits, whereas on field 
'f there is little difference between A and A s . 

The discrepancy between the two fields A s — A appears 
to be caused at least in part by the poorer seeing on parts 
of field 'f. Figure 8 shows that at R = 21.0-23.0 limits, the 
field 'f uj(9) falls at 9 < 3.2 arcsec, indicating that pairs 
are not being resolved, whereas on field 'e' (Figure 7), this 
decrease occurs only at 9 < 2 arcsec. Hence it is field 'e' that 
will provide the more accurate measure of the true excess of 
close pairs. 



5 uj{9) RESULTS 

5.1 Full- field estimates 

Table 2 gives the fitted A and A 3 amplitudes for the two 
mosaiced fields, for galaxies from R = 18.5 to a series of 
faint limits from 7? = 21.0 to R = 23.5. The errors were 
estimated by fitting amplitudes to our measurements of to{9) 
from datasets with one image excluded, and multiplying the 
scatter between these amplitudes by 5.00 for field 'e' and by 
4.25 for field 'f. 

However, the fitted amplitudes are likely to be underes- 
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Figure 7. Observed lu(9) for galaxies on field 'e' brighter than R = 21.0, 22.0, 23.0 and 23.5, as log-linear (left) and log-log (right) plots. 
The dotted lines show the best-fit two-part 0~ - 8 power-laws with the integral constraint offset. 
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Figure 8. Observed for galaxies on field 'f fainter than R = 19.0 and brighter than R = 21.0, 22.0, 23.0 and 23.5, as log-linear 

(left) and log-log (right) plots. The dotted lines show the best-fit two-part 0~ o s power-laws with the integral constraint offset. 
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Table 4. Best-fitting power-laws (5) for io(8) 



R limit 


Field 'e' 


Field T 




8 > 5 arcsec 


8 > 5 arcsec 


21.0 


0.95 ±0.14 


0.84 ±0.13 


21.5 


0.90 ±0.12 


0.88 ±0.20 


22.0 


0.91 ±0.12 


0.87 ±0.18 


22.5 


0.98 ±0.12 


0.87 ±0.18 


23.0 


1.08 ±0.15 


0.88 ±0.18 


23.5 


1.07 ±0.18 


0.90 ±0.20 



5.2 Slope of lo{9) 

Infante and Pritchet (1995) found the slope of u(6) from 
their photographic survey to flatten from 8 = 0.94 ± 0.04 
to <5 = 0.65 ± 0.05 between limits of R = 21 and R = 23, 
and attributed this to a flatter £(r) for later Hubble types. 
However, Couch et al. (1993), using CCD data, found no 
significant change in the u>{9) slope over the same magnitude 
range. We estimate a best-fit slope 5 by least-squares fitting 
l A(9~ s — Cy to u>(9) at 6 > 5 arcsec, for a range of slopes, 
each time recalculating the integral constraint as 



EiV rr ((9) 



(9) 



and finding the values of 8 which minimize the \ 2 of the 
fits. Table 4 gives the best-fit 8, which show no significant 
change over this magnitude range, on either field. If any- 
thing, the power-law is slightly steeper than 8~°' s , although 
this might be due to some of the effect of interactions be- 
tween close pairs of galaxies extending to 9 > 5 arcsec scales. 
The assumption of a 8 — 0.8 slope in fitting an amplitude 
remains consistent with observations within 2a. 



5.3 Individual frame ui(9) 

We also calculate an individual frame u>(9), using the meth- 
ods of Section n except that the 27 CCD frames of field 'e' 
and the 20 of field 'f ' are treated as independent fields, with 
N r = 20000 randoms on each. The uj(9) of the CCD frames 
within each field are averaged to give Lu in (9) for each of the 
two fields. Amplitudes Ai n were then obtained by fitting the 
Uin{9) at separations 5 arcsec < 9 < 21 arcmin with func- 
tions A in (9^°' li — C), where C = 8.28, the larger integral 
constraint factor calculated with equation (8) for a single 
CCD area. 

We find that A in tends to be slightly lower than the 
corresponding A for the full mosaiced fields (Tables 2 and 
3), although the differences are less than ~ 2a, e.g. at the 
R = 23.5 limit, A in = 13.58 ± 1.44 x 10~ 4 for field 'e' and 
A in = 9.33 ± 2.00 x 10~ 4 for field 'f (uncorrected for star- 
contamination, with errors from the scatter between the Ai n 
of the individual CCD frames within each field) 

Figure 9 compares u)i n (9) with the full-field uj(9) at this 
limit. To correct for the difference in their integral con- 
straints, we add to u) in (9) a positive offset of A in (Ci n — 
Cfuii), where d n and Cf r are the integral constraint con- 
stants for single frames and for the full fields respectively, 
giving corrections of 0.00879 for field V and 0.00582 for field 
'f at R — 23.5. The plot also shows the difference of the two 
estimates, Au>(9). 



(a) 



» 1 



-i | i i i i r 
Field 'e' R<23.5 



.full field 

□ individual frame 



(b) 



' i 1 1 1 1 r 

Field T R<23.5 



ill f 



-2 -1 
Log 6 (degrees) 



Figure 9. Full-field w(0) and individual frame w;„(0), the latter 
with offsets added to correct for the difference in integral con- 
straint, for 18.5 < R < 23.5 galaxies on (a) Field 'e' (b) Field 'P. 
The dashed lines show the difference between the two estimates 
Au>(8) = ui(B) — u>i n (0), as a function of log 8. 



Firstly, we note that, at small separations of 9 < 0.01 
deg, the errors on ui(9) derived from the scatter between the 
individual frame estimates are very similar to the error bars 
on the full-field ui(9) estimated by the 'jackknife' method in 
Section 4, as expected if the 'jackknife' method is valid. 

Secondly, Au>(9) may give an indication of the scatter in 
photometric zero-points between the CCD frames forming 
each field. Although some CCD frames were corrected for 
photometric offsets in Section 2.5, there will inevitably be 
some uncorrected scatter a, which will have little effect on 
uJin(9) but increases the frame-to- frame variance in number 
counts by (iVo^ln 10) 2 , where 7 is the number counts slope 
(the PLE model gives 7 = 0.39 at R ~ 23.5) and N the 
mean number of galaxies per frame. This extra variance will 
then increase the full-frame oj(9) at 9 less than the frame 
size by A{a) ~ (a7ln 10) 2 ~ 0.90cr 2 . 

However, our results suggest that, at least on field 'e', 
Au>(9) is not a constant, but contains a ^-dependent com- 
ponent of similar slope to co(9) itself. A likely explanation 
is that the galaxy to(9) amplitude tends to be higher in ar- 
eas of the data with large-scale (similar to or larger than 
the frame size but smaller than the field size) density en- 
hancements, such as rich clusters. In the individual frame 
lo(9) estimates, each frame area is given the same weighting 
whether or not it contains a cluster, whereas the full-field 
estimate would give higher weighting to areas with a high 
galaxy density. The resulting Au)(6) would then contain a 
power- law component in addition to a constant, A (a), from 
the photometric scatter. 

To separate the two components we fit Au)(9) with func- 
tions Aa9~°' s + A(a). At R = 23.5, this gives A A = 1.94 x 
10~ 4 A(a) = 8.26±9.17x 10~ 4 for field 'e', A A = 1.24x 10~ 4 
A(ct) = 1.095 ±0.094 x 10~ 2 for field 'f. On both fields, A A 
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Table 5. Field 'f ui(0) amplitudes at 8 > 5 arcscc (A) and 
2 < 6 < 5 arcscc (A s ) with an estimated correction for the ef- 
fects of photometric scatter, in addition to the correction for star 
contamination. 



R limit 


A (corrected) 


A s (corrected) 


21.0 


30.08 ±3.74 


36.14 ±8.97 


21.5 


20.00 ±4.23 


20.32 ±4.51 


22.0 


18.31 ±3.26 


20.28 ±4.20 


22.5 


15.73 ±3.83 


23.70 ± 3.44 


23.0 


16.49 ±3.89 


25.33 ±4.30 


23.5 


18.59 ±4.39 


27.63 ±4.35 



is about an order of magnitude less than the u>(8) amplitude, 
and therefore comparable to the ui(8) errors and cannot be 
regarded as an accurate measurement. However, a correla- 
tion between ui(6) on 9 < 0.1 deg scales and galaxy density 
might more effectively be quantified using the moments of 
the galaxy counts-in-cells. As uj{9) is related to the variance 
of the galaxy counts-in-cells above the Poisson expectation, 
the higher variance in areas of the data with higher counts 
will will result in the skewness and kurtosis also exceeding 
the Poisson expectation - we investigate these higher mo- 
ments in Section 8. 

We estimate the photometric scatter a as y 7 , giving 
a = 0.030±^3o mag for field 'e' and a = 0.110 ±6.005 mag 
for field 'P. To estimate the likely effect on our results, we 
fit amplitudes to our full-field uj(9) with offsets subtracted 
at 9 < 0.3 deg, of (§§) x 8.26 x 10~ 4 for field 'e' and (±§) x 
1.095 x 10 -2 for field 'f (the factors in brackets allow for the 
partial cancelling out of the negative offset by the reduction 
in integral constraint, smaller than the offset by the ratio 
of the frame to the full-field area) ; these offsets reduced the 
u>{9) amplitude at 9 > 5 arcsec by 1.1 per cent on field 'e' 
and 21.2 per cent on field 'P. 

Hence the effects of photometric scatter on ui (9) appear 
to be quite insignificant on field 'e', but on field 'P are suf- 
ficient to cause a > la overestimation of the amplitude at 
the faint limits of the data, so some sort of correction should 
be applied. We fit amplitudes A and A s to the field 'P lo(9) 
with the negative offset described above, at all magnitude 
limits, and apply the same star-contamination corrections 
as previously, giving a new set of field 'P lu(9) amplitudes 
(Table n) with corrections for both star-contamination and 
photometric scatter. 

This correction reduces the large scale A much more 
than A s , so that some difference in these amplitudes is now 
seen at R > 22.5, but only at the R = 23.5 limit does A s — A 
become similar to that on field 'e'. This suggetsts that pho- 
tometric errors can account for part of the difference in the 
uncorrected A 3 — A of the two fields, but the effects of poorer 
seeing on field 'P remain of equal or greater importance. For 
the remainder of the paper, we assume these new amplitudes 
for field 'P, but note that such corrections for the effects of 
photometric scatter can only be estimates and we must re- 
gard the field 'P results as less reliable than those from field 
'e\ 

5.4 lo(9) scaling with magnitude limit 

Figure 10 shows the scaling of the u)(9) amplitude with R 
magnitude limit, for this data and several other 7?-band 



surveys, plotted as the amplitude at one degree when a 
(9~ 0,8 power-law is fitted. The magnitude limits of other sur- 
veys were converted into approximate equivalents in our R 
band. For the four red-band CCD surveys, we assume R = 
i?(Roche et al. 1993) -0.02, R = R( Roche et al. 1996) ±0.02, 
R = R (Woods and Fahlman) and R — r(Brainerd et 
al.) — 0.55, and for the Infante and Pritchet (1995) F-band 
photographic survey, R = F~ 0.14 from the correction given 
by Metcalfe et al. (1991). Couch et al. (1993) used a broad 
'VR' passband, and Villumsen et al. (1997) give u>(9) am- 
plitudes from the Hubble Deep Field as a function of Vme 
limit. For these bluer passbands, we assume R = VR — 0.35 
and R = Veo6 — 0.39, derived using the modelled spectrum 
of an evolving Sbc galaxy (Section 6.1) at z ~ 0.5. 

There is some scatter between the observations, with 
our u){9) amplitudes, especially those from field 'e', sup- 
porting the relatively high normalization of the Infante and 
Pritchet (1995) and the Woods and Fahlman (1997) u(6) 
scaling. The results from this data alone might suggest a 
levelling-out of the to(9) scaling at R ~ 22.5-23.5, but a com- 
parison with Brainerd et al. (1995), Metcalfe et al. (1995) 
and Villumsen et al (1997) indicates that, at least at these 
A ~ 0.4-0.8/im passbands, the to(9) amplitude falls much 
further beyond our survey limit. The interpretation of these 
results is discussed further in Sections 6.3 and 9. 



6 MODELS OF lo(9) 

6.1 Modelling of Galaxy Evolution 

In order to interpret the faint galaxy ui(9), we need a model 
which generates a redshift distribution N(z) for each type 
of galaxy, which is as consistent as possible with observa- 
tions. In this paper the results are compared with a Pure 
Luminosity Evolution (PLE) model, i.e. a model in which 
the characteristic luminosity of the galaxy luminosity func- 
tion, L* , evolves with redshift but the normalization <f>* and 
slope a are assumed constant. 

In the 'size and luminosity evolution' model of Roche 
et al. (1998a), elliptical galaxies were assumed to form in 
a short starburst at high redshift but spirals form by a 
more gradual inside-outwards process. Both elliptical and 
spiral galaxies undergo an evolutionary brightening out to 
z ~ 1-2. The PLE model of this paper is an updated ver- 
sion of this in which (i) the spectral evolution is now mod- 
elled using new 'bc96' stellar evolutionary models (e.g. Char- 
lot, Worthey and Bressan 1996), with a Salpeter IMF, solar 
metallicity for E/SO/Sab/Sbc galaxies and 0.4 solar metal- 
licity for Scd/Sdm and Starburst galaxies, with the same 
star-formation histories as previously (ii) dust extinction in 
the galaxies is now modelled as an optical depth 

a 10 10 dM 
Td ust (A) = Pdust X — rr™^) (10) 

Mgai at 

where ^ is the SFR in M© yr -1 , M ga i the total galaxy 
mass (stars and gas), and tm(A) the tabulated 'Rv=3.1' 
Galactic dust extinction law of Mathis (1990), normal- 
ized to unity at 4500A, and Pdust = 1-0 for E/S0 galax- 
ies, 0.4 for Sab spirals 0.3 for Sbc spirals and 0.12 for 
Scd/Sdm/Starburst galaxies. 
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Figure 10. The u>{6) amplitude at 9 > 5 arcsec for the galaxies in our two fields, against R magnitude limit, compared with ui(9) 
amplitudes from other faint galaxy surveys and five models described in Section 6; model A, non-evolving, e = (dotted); model B, PLE, 
e = (solid); model C, PLE, e = 1.2 clustering evolution (short-dashed); model D, PLE e = —1.2 comoving clustering (long-dashed); 
model E, PLE, e = 0, with a L -0 ' 25 decrease in the strength of clustering with luminosity (dot-dash). 
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Table 6. Luminosity function normalizations <E>* (in galaxies 
Mpc -3 for Ho = 50 km s _1 Mpc _1 , slopes a and characteristic 
_R-band absolute magnitudes M* (for Ho = 50 km s _1 Mpc _1 ) 



Galaxy type 


$* 


a 


M R 


E 


0.48 x 10~ 4 


+0.54 


-22.14 


SO 


8.6125 x 10~ 4 


-0.12 


-22.09 


Sab 


10.0675 x 10~ 4 


-0.32 


-21.76 


Sbc 


9.1625 x 10~ 4 


-0.71 


-21.71 


Scd 


2.125 x 10~ 4 


-1.45 


-21.89 


Sdm 


1.10625 x 10~ 4 


-1.89 


-21.87 



The galaxy luminosity functions in the model are those 
derived by Bromley et al. (1998) from the red-band Las 
Campanas redshift survey, divided into six spectral types. 
The spectral 'clans' 1, 2, 3, 4 and 5 are assumed to corre- 
spond to our PLE models for E, SO, Sab, Sbc and Scd types 
respectively, while half of the galaxies in the bluest class 
(clan 6) are described with our evolving Scd model (these 
could represent Sdm/Irr galaxies in a non-starburst phase) 
and the other half with a 1 Gyr age starburst k-correction. 
We adopt the suggested 25 per cent for the fraction of light 
missed in this survey's isophotal magnitudes, asnd hence 
apply a correction of -0.31 mag to the Mr, and also incor- 
porate the estimated corrections to the luminosity function 
slopes for incompleteness due to surface brightness biases. 

Table 6 lists the parameters of the luminosity func- 
tions, which are Schechter functions of the form <j>{L) = 
<&*(L/L*y a , used in our model. 

For comparison, we also consider a non-evolving model, 
with the same z — luminosity functions, and k-corrections 
derived by redshifting without evolution the galaxy spectra 
of the PLE models at z = 0. 

The L* evolution in the PLE model increases the pre- 
dicted galaxy number counts sufficiently to fit the obser- 
vations(Figure 5). Figure 11 shows the galaxy redshift dis- 
tributions N(z) predicted by both models at R = 21-23.5 
limits. In the PLE model, a 'tail' of 1 < z < 2 galaxies 
appears between R = 22 and R = 23, whereas with no L* 
evolution very few z > 1 galaxies would be seen even at the 
survey limit of R = 23.5. Deep spectroscopy (Cowie et al. 
1996) reveals that 1 < z < 1.8 galaxies do start to appear 
in substantial numbers at B > 23.2 (R <; 22.7), in general 
agreement with this PLE model. The PLE model is also rea- 
sonably consistent with the numbers and surface brightness 
of the z > 2.5 'Lyman break' galaxies on the Hubble Deep 
Field, although these more distant galaxies are not expected 
to be visible within this survey. 

Figure 12a shows the mean redshift of both models as 
a function of R magnitude limit. At the survey limit of 
R = 23.5, the PLE model predicts z mean — 0.685, which in- 
creases on going further faintward to z mean ~ 1.4 at R ~ 27. 
As our PLE model seems to fit the galaxy counts and N(z) 
reasonably well, we now use it in interpreting the w(0) am- 
plitudes. 

6.2 From N(z) to w(0) 

The u>(6) amplitude is predicted using the Limber's formula 
integration of £(r, z) over the modelled N(z). We assume the 
form given in the Introduction, 




Redshift 

Figure 11. Galaxy redshift distributions predicted for all galaxies 
to limits of (from lowest to highest), R = 21.0, 22.0, 23.0 and 23.5, 
for the PLE (solid) and no-evolution (dotted) models. 

£(r, Z ) = (r/r )^(l + z)- (3+E) (11) 

where ro normalizes the strength of clustering at z = 0, 
7 is the slope and e parameterizes the clustering evolution 
relative to the e = stable clustering model. 

The uj{6) amplitude will generally decrease on going 
faintward, primarily as a result of N(z) becoming more ex- 
tended, but there may be additional effects due to the evo- 
lution of clustering and to changes in the proportions of 
different types of galaxy. The model of Roche et al. (1996) 
assumed 7 = 1.8 for all galaxy types, with ro = 5.9 h^ 1 
Mpc for E/S0 galaxies and ro = 4.4 h^ 1 Mpc for later types 
(from Loveday et al. 1995, hereafter LMEP) and also dwarf 
galaxies (those with 2 = blue-band absolute magnitudes 
Mb > —20.5) half as clustered as those of higher luminosity 
(again from LMEP). 

With stable clustering and L* evolution, this model pre- 
dicted a u)(9) scaling fairly consistent with observations in 
the blue-band (Roche et al. 1996), but it underpredicted 
the u>{6) amplitudes from Jf-band surveys at K ~ 19.5- 
21.5 limits (Roche et al. 1998b; Carlberg et al. 1997). The 
results from the Jf-band surveys, in which E/S0 galaxies 
will be more prominent, were fitted better by assuming even 
stronger clustering for early-type galaxies (e.g. as observed 
by Guzzo et al. 1997). 

The ro used by Roche et al. (1996) were derived by 
LMEP from a cross-correlation of each type of galaxy with 
all APM galaxies, a method which would tend to underes- 
timate the differences between the £(r) of the galaxy types. 
LMEP also derive £(r) using a different 'inversion' method, 
which gave r = 7.76±0.15 h' 1 Mpc and 7 = 1.87±0.07 for 
E/S0 (of all luminosities) and r = 4.49±0.13 h' 1 Mpc and 
7 = 1.72±0.05 for late-types. These estimates would be more 
consistent with Guzzo et al. (1997) and the results from the 
7f -band. In this paper we model £ (r) with, for simplicity, a 
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Figure 12. (a) Mean redshift as a function of R magnitude 
limit for PLE (solid) and non-evolving (dotted) models, (b) 
The fractions of E/SO galaxies and dwarf galaxies (those with 
Mb > —20.5 at z = 0) as a function of R magnitude limit for 
the PLE model, (c) The mean Mr at z = as a function of 
magnitude limit in the PLE model. 



single slope 7 = 1.8, and take the normalization from the 
LMEP 'inversion' estimates at r = l/i -1 Mpc. This gives 
for E/SO galaxies (ro) 18 = 7.76 1 ' 87 and thus r = 8.4 h' 1 
Mpc, and for later types (ro) 1 ' 8 = 4.49 1 ' 72 giving r = 4.2 
h- 1 Mpc. 

In the Limber's formula integration, we normalize our 
models to the LMEP spiral ro and then apply a weighting 
term to the galaxy clustering at each redshift derived from 
the predicted fraction of ellipticals (fe)- With our assumed 
ro, ^elliptical — 3-5£s P irai, and if the cross-correlation be- 
tween early and late type is assumed to be the geometric 
mean of their autocorrelations so that Across — l-87£ sp j ro i, 
the adopted weighting factor is (1 — / e ) 2 + 3.5/1 + 2 x 
1.87/ e (l - fe) - or equivalently (1 + 0.87/ e ) 2 . 

The luminosity dependence of clustering is modelled by 
two methods. The first is as in Roche et al. (1996) - accord- 
ing to LMEP, lower luminosity (Mb > —20.5) galaxies are 
less clustered than more luminous galaxies of the same Hub- 
ble type by a factor of 2 at r ~ 1 h~ Mpc, so, at each red- 
shift, the models calculate fd, the fraction of galaxies which 
would have Mb > —20.5 at 2 = 0, and apply a luminosity 
weighting, normalized to be unity for the dwarf fraction in 
the LMEP dataset, of 1.2(l-/ d ) 2 +0.6/| + 2 x 0.85(l-/ d )/ d . 

Figure 12b shows f e and fd as a function of R limit in 
the PLE model - f e slowly decreases (from 0.28 to 0.18 be- 
tween R < 21 and R < 23.5) while fd increases. Both these 
trends would slightly steepen the scaling of ui(9). Dwarf 
galaxies (almost all of late type, at least in the field environ- 
ment) will dominate the sample at R <i 26, as L* galaxies are 
seen to their highest redshifts at these magnitudes and going 
fainter simply results in looking further down the faint-end 



of the galaxy luminosity function for the same volume of 
space. 

Using the luminosity evolution and ro described above, 
we predict the scaling of u>(0) with R magnitude limit for 
four models 

A: N(z) from non-evolving model, e — 
B: N(z) from PLE model, e = 

C: N(z) from PLE model, clustering evolution of e = 1.2. 
D: N(z) from PLE model, comoving clustering of e = —1.2. 

We also consider a fifth model E, which differs from 
model B only in the luminosity dependence of £(r). Here 
the strength of clustering varies continuously with the un- 
evolved luminosity (i.e. very approximately the mass) of the 
galaxies, as £(r) oc L ' 25 . The luminosity weighting term in 
Limber's formula is then io- 01 « M R>+ 22 -°) ; where (Mr) is 
the mean unevolved (i.e. at z = 0) .R-band absolute mag- 
nitude of all galaxies at a given redshift, and the addition 
of 22.0 normalizes the model to the same u>(6) amplitude as 
model B at R ~ 22.0. Model E would be consistent with 
the dependence of £(r) on absolute magnitude as observed 
(LMEP, Guzzo et al. 1997) for relatively luminous galaxies, 
but includes the additional assumption that clustering con- 
tinues to decrease below Mb = —20, becoming very weak 
for the faintest dwarfs. Models B and E would be indistin- 
guishable at the limits of our survey but differ greatly at 
7? <; 26 where the mean unevolved luminosity (Figure 12c), 
which decreases only slowly to R ~ 25.5, falls more steeply 
as the sample becomes dominated by dwarf galaxies. 



6.3 Comparison with Observations 

Figure 10 shows the five models against the u>(0) observa- 
tions. Model A gives co(0) amplitudes well above most of 
the data and can be rejected (redshift surveys have already 
excluded models with no L* evolution). To R ~ 25, PLE 
model B passes through the middle of the rather scattered 
data points and is the most consistent with our results; how- 
ever it overpredicts the very weak clustering of the faintest 
Hubble Deep Field galaxies. 

Model C, with clustering evolution of e = 1.2, passes 
through only the lowest of the data points at R < 25. At 
R = 21-23.5, e = 1.2 clustering evolution reduces ui(6) 
0.134-0.191 dex below the model B prediction. It signifi- 
cantly underpredicts our amplitudes for field 'e' at all limits 
and for field 'f at R = 23-23.5, and appears to be rejected 
by > 3a. The comoving model D gives a very similar lo(9) 
scaling to the non-evolving model A, overpredicting the faint 
galaxy clustering, and can be rejected. 

If we assume the N(z) and ro in our PLE model to be 
correct and treat e as a free parameter, a comparison of the 
model at R = 21-23.5 with the combined field 'e' and field 'f ' 
results minimizes X 2 fore = 0.02l^ 8 .The same comparison 
for the results from the more photometrically reliable and 
less star-contaminated field 'e' only gives a consistent e = 
-0.29±° l\, and rejects e > 0.77 and e < -1.03 by > 3a. 

To R ~ 24, the model E uj(6) does not differ signifi- 
cantly from model B, so it is similarly consistent with our 
results, but at R > 26 the decrease in ro with luminosity 
causes the u>(9) amplitude to continue falling after z mean 
has reached its maximum, in contrast to models A-D where 
the increasing number of low redshift dwarfs causes lo(9) 
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to start increasing again. Model E, in contrast to the other 
four, fits the Hubble Deep Field clustering very well. 



7 THE SMALL-SCALE w(0) AND MERGER 
RATE EVOLUTION 

Infante et al. (1996) found that the u>(6) of a sample of 
R < 21.5 galaxies, remained a S ~ 0.8 power-law at 9 < 6 
arcsec but with an amplitude a factor 1.8 higher than at 
9 > 6 arcsec. We investigated this effect in our data (Section 
4) by fitting two separate ui(9) amplitudes A and A s at 9 > 5 
arcsec and 2 < 9 < 5 arcsec respectively, and found that on 
field 'e' A s was significantly higher than A, at all magnitude 
limits. The results from field 'f ' are not used in the following 
analysis, as although the corrected amplitudes for this field 
did show some excess of A s over A at the faintest limits 
(Table n), any estimate of A s — A will be very dependent on 
our uncertain correction for photometric scatter, and even 
after this correction A s would be, at least at R < 23, affected 
by the poorer seeing. 

At the mean redshift of the Infante et al. (1996) sample, 
Zmean — 0.35, 6 < 6 arcsec corresponds to a projected phys- 
ical separation ^ 20 /i -1 kpc. Even locally, there is an ex- 
cess of <, 20 h- 1 kpc pairs above the number expected from 
normal galaxy clustering, made up of interacting galaxies in- 
cluding those in the early stages of merging. In a magnitude- 
limited (B < 14.5) = 0.007 subsample of the UGC 
catalog, 4.6 ±0.4 per cent of galaxies are in physical pairs of 
projected separation < 19 ftT 1 kpc (Carlberg et al. 1994). 
We derive a similar 'close pair fraction', f pa irs, for our field 
'e', here defining fpairs as the fraction of the galaxies which 
are in 6 < 5 arcsec pairs, above the expectation from the 
lo(8) measured at 6 > 5 arcsec. 

Firstly, we assume that ui(9) oc 0~ 8 at all 9 < 5 arcsec, 
so that we can estimate the number of 9 < 2 arcsec pairs 
from an inwards extrapolation of this power-law with its 
2 < 9 < 5 arcsec amplitude. In an unbounded area, with a 
surface density of galaxies p ga i, the number of 9 < /3 pairs 
per galaxy above the random expectation will be 



N pair 

N, 



gal 



PgalA 



(12) 



where A is the u>(9) amplitude. To find the number of pairs 
above that expected from the lo(9) seen at larger scales, A is 
replaced by the difference in clustering amplitudes A s — A, 
so that 



N„ 



Pgal(A s - A) 



Jo 



Ngal 

which for f3 = 5 arcsec gives 

N pair 



2ir99-° s d9 



Ngal 



0.00196p ga i(As - A) 



(13) 



(14) 



In a finite field N pa irs will be reduced as some galaxies will 
lie less than 5 arcsec from the field edges or the holed areas, 
but this effect is minor for the large area of our fields - a 
numerical integration over the field 'e' area gives 



N pair 



gal 



= 0.001938/95^(^4.3 — A) 



(15) 



Table 7. Percentage of galaxies in excess 6 < 5 arcsec pairs on 
field 'c' 



R limit 


fpairs (percentage) 


21.0 


4.03 ± 1.85 


21.5 


5.14 ±2.08 


22.0 


8.32 ±2.60 


22.5 


12.33 ±3.36 


23.0 


12.88 ±4.05 


23.5 


11.12 ±4.86 



As a pair contains two galaxies, f pa ir = 2 . We esti- 

mate fpairs at each magnitude limit using the field 'e' A 
and A s corrected for star contamination (Table 3), and also 
apply a star contamination correction to the observed p ga i, 



fpa t rs = 0.003876(1 - fstar)Pgal(As - A) 



(16) 



which is given in Table 7, with errors from adding the errors 
on A and A s in quadrature. The excess of close pairs above 
the u>(9) measured at larger scales is 3.7cr significance at 
R < 22.5 and > 3cr at R < 22 and R < 23. 

Infante et al (1996) found 1317 pairs with 2 < 9 < 6 arc- 
sec in a sample of 16749 galaxies with R < 21.5, compared 
to 477 expected by chance or 842 expected from the cluster- 
ing at larger scales. The proportion of galaxies in 2 < 9 < 6 
pairs above the large-scale uj{9) was then 2 x 13 |g^| 42 = 5.67 
per cent. For a 9~°' s power-law the number of 9 < 5 arcsec 
pairs will be 1.097 times the number at 2 < 9 < 6 arcsec, and 
the significance of the pairs excess is given as ~ 5a, so the 
Infante et al. (1996) result corresponds to f pa i r = 6.22±1.24 
per cent. 

As the timescale between the approach of two galax- 
ies within ~ 20 h" 1 kpc and their subsequent merging or 
separation, ~ 0.75 Gyr (e.g. Mihos and Hernquist 1996), is 
short compared to the Hubble time, the merger/interaction 
rate R m oc f pa i r approximately. The merger rate evolution 
is modelled as R m oc -R m o(l + z) m and normalized to a 4.6 
per cent fraction of local galaxies in pairs of projected sepa- 
ration < 19 h^ 1 kpc (Carlberg et al. 1994). Pairs within this 
separation will have an angular separation 9 < 5 arcsec at 
angular diameter distance dA > 784 h^ 1 Mpc (z > 0.455 in 
our chosen cosmology). At dA < 784 hT Mpc, if we again 
assume u>(9) oc 9~ ' 8 for the close-pair galaxies, the fraction 
with 9 < 5 arcsec will be ( 7 ^1 ) 1 ) 1 - 2 . 

For each magnitude limit we then model f pa ir by sum- 
ming 

fpair = fpair(z < 0.455) ± f pal r(z > 0.455) (17) 

where 



f P air(z < 0.455) = 0.046 



Jo°- 455 (7it^r) 1 - 2 ^)( 1 + ^ m ^ 
C 455 N(z)dz 



fpair(z > 0.455) 



f. r N(z)(l + z) m dz 

QM6 J0455 >_ 



/(MBB 

over the PLE model N(z). 

Infante et al. (1996) estimated m = 2.2 ± 0.5 from the 
increase in the fraction of galaxies in close pairs from z = 
to z mean = 0.35. Our data suggest that the upward trend in 
fpair continues to fainter magnitudes, with the slight drop 
from R = 23 to R = 23.5 probably being due to less effective 
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Figure 13. Fraction of galaxies in 9 < 5 arcsec pairs in excess 
of the number expected by chance and the clustering observed 
at 6 > 5 arcsec, as a function of R magnitude limit, from our 
field 'e' and from Infante et al. (1996). The dotted line shows a 
model in which the local fraction of galaxies in close pairs remains 
constant, the dashed lines are models in which the merger rate 
evolves as (1 + z) m , where m = 1, 2, 3, 4. 



splitting of close pairs on approaching the limits of the data 
(fpair decreases steeply beyond the R = 23.5 limit). 

Figure 13 shows f pa ir from our field 'e' and from In- 
fante et al. (1996), compared with the model (equation 18) 
for m — 0,1,2,3 and 4. The m = 2 model is the most 
consistent with the data; both a non-evolving merging rate 
(m = 0) and very rapid evolution of m = 4 appear to be 
rejected. Our R = 21.0 point lies significantly below the 
m — 2 model, but at the brighter limits there might be 
some underestimation of f pa ir if the larger angular sizes of 
brighter galaxies prevents the resolution of some 2 < 6 < 5 
arcsec pairs. A \ 2 test of our model against the plotted data 
(excluding the R = 21.0 point) gives a best-fitting merger 
rate evolution of m = 2.01^g'gg, with evolution of m > 3.25 
rejected by > 3a. This is discussed further in Section 9.2. 



8 HIGHER-ORDER CORRELATION 
FUNCTIONS 

The large size of our galaxy sample enables us to investigate 
their higher order correlation functions, which may discrim- 
inate between models of galaxy clustering more effectively 
than ui(6) alone. The higher order correlations have previ- 
ously been studied at brighter magnitudes of 17 < B < 20, 
for 1.6 x 10 6 galaxies in the APM survey (Gaztanaga 1994) 
and 2.9 x 10 5 galaxies Edinburgh-Durham Southern Galaxy 
Catalog Survey (EDSGCS) (Szapudi, Meiksin and Nichol 
1996), in the form of hierarchical moments. Before discussing 
these results further, we define the hierarchical moments 



and describe their estimation from a simple counts-in-cells 
method. 



8.1 Estimation of Hierarchical Moments 

Consider a survey area divided into m cells of area fl, with 
Ni galaxies in cell i and a mean number of galaxies per cell 
N. The moments of the distribution of counts-in-cells will 
be 



m 



(18) 



and the second moment (variance) of the counts-in-cells will 
be related to u>(6) as 



/j 2 = N + N 2 



uj(e)dQ 1 dn 2 



(19) 



where the double integral is over the area of one cell. An 
area-averaged u>(6), C02, can be defined as 



Ul2 



Lu(e)dn 1 dn 2 



(20) 



and estimated from the measured second moment by sub- 
tracting the expectation for a Poission distribution (i.e. 
fJ-2 = N) and dividing by iV 2 . 

U2=^2-N) (21) 

Similarly, the area-averaged angular three-point correlation 
function 0)3 can be calculated from the excess of the third 
moment (skewness) relative to the Poission expectation 
(Baugh, Gaztanaga and Efstathiou 1995) 



^3 = ^(M3-3/i 2 + 2iV) 



(22) 



and the area-averaged angular four-point correlation func- 
tion W4 can be calculated from the excess of the fourth mo- 
ment (kurtosis) relative to the Poission expectation 



UI4 



3fxl + 11/42 +6iV) 



(23) 



It is useful to express higher order correlation functions in 
the form of hierarchical amplitudes sj, defined as 



sj 



OJJ 



(24) 



We divide fields 'e' and 'f ' into square cells of side 9, 
and count the number of 18.5 < R < 23.5 galaxies in each 
cell. Many of the cells will overlap with field edges or holed 
areas. If the missing proportion of the cell area f m is less 
than 0.3, the galaxy count is corrected by multiplying by 
(1 — /m) _1 , while cells with f m > 0.3 were excluded from 
the analysis (as in Gaztanaga 1994). Combining fields 'e' 
and 'P, we calculate the variance, skewness and kurtosis of 
the counts-in-cells, for cell size 8 ranging from 0.01 to 0.1 
degrees. 

As for u)(6), u>2(8) is corrected for star-contamination 
by multiplying by (1 — /s) -2 , where the estimated star- 
contamination f s = 0.15866 ± 0.00634 for fields 'e' and 'f 
combined. Error bars were calculated by randomly divid- 
ing the cells into 10 subsamples, calculating 012(0) for each 
subsample the scatter between them by ^/10, finally adding 
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Figure 14. The area-averaged angular correlation function u>2(8), 
from the counts-in-cclls analysis of both fields combined, with the 
best-fitting function A(6»-°' 8 - 0.9606). 



in quadrature the errors from the uncertainty in the star- 
contamination correction. 

As for ui(6), an amplitude is estimated by fitting with a 
function A(8~ 0,s — C), where the integral constraint factor 
C is derived from those of the two independent fields, C e 
and Cf (see Section n) and their areas S7 e and Q,f, as 



&iC e + SifCf 



(25) 



giving C — 0.9606. Figure 13 shows ui2(8) with the best- 
fitting function A(6-°- 8 - 0.9606), with A = 56.8±1.3 x 10~ 4 
at a cell size one degree. As the integral constraint factor 
C — 2.24 for a 1 x 1 deg cell, this corresponds to a ui(8) 
amplitude of = 25.4 ± 0.6 x 10~ 4 , slightly higher than 
but reasonably consistent with the amplitudes calculated 
from the positions of individual galaxies. 

To correct the hierarchical moments for the integral 
constraint in u)2(6), estimated as AC = 0.00387 from the fit- 
ted 012(6 amplitude before correction for star-contamination, 
83(6) and sa(6) are calculated as 



sj(0) 



(uj 2 (8) +0.00387)- 



(26) 



No corrections are applied for the integral constraints in 
LJ3,{9) and Qa(6), but these will be much smaller than for 
0)2(8), due to the much steeper slopes of the higher-order 
correlation functions. 

Star contamination will reduce 013(6) by a factor (1 — 
/ 3 )~ 3 but reduces the denominator u)2(6) 2 by (1 — /s)~ 4 , 
so 33(6) is actually overestimated, and can be corrected for 
star-contamination by multiplying by (1 — fsr)- Similarly, 
star contamination will reduce uj±(6) by (1 — / s )~ 4 but re- 
duces ui 2 (6) 3 by (1 - / s )" 6 , so sa(6) is corrected by multi- 
plying by (1 - f s f. 



By a similar argument, frame-to-frame photometric 
scatter increases u) 2 (6) over this range of cell sizes, and so is 
expected to decrease 83(6) and 84(8), with the effect on 83(6) 
being about half that for 10(6). The photometric scatter on 
field 'f was estimated to reduce the field's u)(6) amplitude, 
fitted over a similar range of 6, by 21.2 per cent (Section 
5.3), so for both fields combined its effect on 83(6) should 
be < 5 per cent, much less than la. However, it is possible 
that photometric or other differences between the two fields 
could also affect results for the combined fields. In order to 
check the validity of using the combined fields, we compare 
the results with 83(6) and 34(6) derived by performing the 
same counts-in-cells analysis for the cells in field 'e' only, ap- 
plying corrections for star-contamination and integral con- 
straint appropriate for this field alone. 



8.2 Results for s 3 (6) and s 4 (0) 

Figures 14 and 15 show 83(6) and 84,(6), corrected for both 
star-contamination and the integral constraint. The s 3 (6) 
and S4(6) of the field 'e' data only are very similar in both 
slope and normalization to those from the full dataset, sup- 
porting the accuracy of the combined field result. There is 
a highly significant (~ 8a) detection of a positive signal in 
83(d), and it is clear that 83(8) is not constant with increas- 
ing 8 but decreases quite steeply over this range of cell size - 
the best-fitting 8~ s power-law is of slope S = 0.413 ± 0.083. 
The statistics are poorer for S4(8) but the detection is still 
4-5cr and the best-fit slope of S = 0.614 ± 0.134 is even 
steeper. 

We compare our results with the 83(6) and S4,(8) of 
17 < B < 20 galaxies in the APM (Gaztanaga 1994) 
and EDSGCS (Szapudi et al. 1996) surveys. The hierar- 
chical moments are projections of the corresponding three- 
dimensional Sj(r), defined as the ratios of the three- 
dimensional J-point correlation functions to £(r) J ~ 1 . Hierar- 
chical moments are well-suited for a comparison of deep and 
shallower surveys, as (i) they vary much less steeply with 
cell size than the uij, and (ii) the projection factors from 
the three-dimensional S,j(r) to the two-dimensional sj(6) 
are almost constant with survey depth. 

As the S3 and S4 projection factors modelled by Sza- 
pudi et al. (1995) remain constant within 4 per cent for 
3 mag shifts in survey limit, we neglect any change in 
sj(6)/Sj(r) between the 17 < B < 20 and 18.5 < R < 23.5 
surveys to be negligable and plot the sj(6) from the shal- 
lower surveys without any vertical shift. However, we must 
take into account the difference in the mean proper distance 
r corresponding to a given cell size 8 in the deep and shal- 
low data. Our PLE model gives a mean angular diameter 
distance of 335/i Mpc for 17 < B < 20 galaxies (with 
Zmean = 0.14) and 860ft _1 Mpc for 18.5 < R < 23.5 galax- 
ies (with z mean = 0.69), so we plot the APM and EDSGCS 
sj(6) shifted horizontally by A(log6>) = log§§§ = -0.41. 
For all data on Figures 15 and 16, log 8 — — 1 will then 
correspond to, on average, r ~ 1.50/i _1 Mpc. 

Szapudi and Gaztanaga (1998) discuss the differences 
between the APM and EDSGCS Sj(8). These are not cor- 
rected for integral constraint, resulting in spurious upturns 
at 8 > 3 deg for the APM survey and 8 > 1.5 deg for the 
smaller-area EDSGCS survey, which account for the differ- 
ence in these two sj(8) at the largest 8 plotted here. How- 
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-1.5 -1 

Log 9 (degrees) 



Figure 15. The hierarchical moment 53(f) of our 18.5 < R < 
23.5 galaxies, for both fields combined and for field 'e' only (offset 
slightly for clarity), compared with the S3 (9) of 17 < B < 20 
galaxies in the APM (Gaztanaga 1994) and EDSCGS (Szapudi 
et al. 1996) surveys (shifted by A(log#) = —0.41 to take into 
account the difference in survey depth). The dashed line shows 
the power- law (6»-°- 413 ) best-fitting our results for the combined 
fields. 



.This work 18.5<R<23.5 
□ 18.5<R<23.5 field V only 

oEDSGCS 17.0<B<20.0 (rescaled in 8) 
a APM 17.0<B<20.0 (rescaled in 8) 



-1.5 -1 

log 8 (degrees) 



Figure 16. The hierarchical moment sn{9) of our 18.5 < R < 
23.5 galaxies, for both fields combined and for field 'e' only (offset 
slightly for clarity), compared with the s^{9) of 17 < B < 20 
galaxies in the APM (Gaztanaga 1994) and EDSCGS (Szapudi 
ct al. 1996) surveys (shifted by A(log0)) = -0.41 to take into 
account the difference in survey depth). The dashed line shows 
the power-law (9 -0,614 ) best-fitting our results for the combined 
fields. 



ever, at 6 < 0.125 deg, here plotted as log(#) < —1.3 and 
corresponding to r < 0.75ft -1 Mpc, the APM and EDSGCS 
sj(9) are genuinely inconsistent, diverging greatly. 

Our results, for both S3(8) and 84(6), are clearly much 
more consistent with the steeper slopes of the EDSGCS mo- 
ments over the same range of physical scales (0.15 < r < 
1.5ft -1 Mpc). There is some indication that our sj(9) are 
lower in amplitude than those from the EDSGCS, but the 
reduction is relatively small - at \og(6) ~ 1.6, the power-laws 
best-fitting our results fall below the EDSGCS moments by 
only 24 ± 13 per cent for s a (6 ) ) and 39 ± 18 per cent for S4,(0)- 
We discuss the implications of the sj(6) slope and amplitude 
in Section 9.3. 



9 DISCUSSION 

9.1 Clustering Evolution 

Using a large sample of ~ 70000 galaxies on two fields, 
we have measured the angular correlation function to a 
magnitude limit of R = 23.5. After corrections for star- 
contamination, our results support the relatively high nor- 
malization of lo(9) from the Infante and Pritchet (1995) and 
Woods and Fahlman (1997) surveys, which used very high 
Galactic latitude fields with minimal star-contamination. 

We can also compare our uj{9) amplitudes with those 
from surveys in other passbands at limits with a similar 
galaxy surface density. Our ui(9) amplitudes at R ~ 22.5- 
23.5 are similar to those measured at / ~ 22-24 (Neuschae- 
fer et al. 1995; Brainerd and Smail 1998) and K ~ 19.5-21.5 
(Carlberg et al. 1997; Roche et al. 1998b) but higher than at 
limits of B ~ 24.5 (Roche et al. 1993) and V ~ 24.5 (Woods 
and Fahlman 1997). This would be expected from the higher 
lo(9) amplitude of redder galaxies when the sample is di- 
vided by V - I (Neuschaefer et al. 1995) B — R (Roche et 
al. 1996) or U -K (Carlberg et al. 1997). Neuschaefer et al. 
(1997) also find the co(8) amplitude of bulge-profile galaxies 
at I < 23 to be 2-4 times higher than that of disk galax- 
ies, suggesting that the stronger clustering of earlier Hubble 
types (LMEP, Guzzo et al. 1997) is maintained to at least 
z ~ 0.7. This might explain why the increase in u(6) to- 
wards longer survey wavelengths appears to occur primarily 
between the V and R passbands, as these would lie on either 
side of the 4000 A break for galaxies near the peak of N(z) 
at these limits, and the size of this break correlates strongly 
with Hubble type. 

Our results are consistent with the strength of galaxy 
clustering measured locally (LMEP), combined with galaxy 
L* evolution consistent with current redshift surveys, and 
galaxy clustering approximately stable (e = 0) in proper co- 
ordinates. Strong clustering evolution of e = 1.2, the linear 
model evolution of e = 0.8, and comoving clustering (e = 
— 1.2) are all disfavoured by ~ 3a. The clustering simulations 
of Colin et al. (1997), for the distribution of mass, predict 
e = 1.08 ± 0.09 with Q = 1 and e = 0.18 ± 0.12 with fi = 
0.2. Hence, on the face of it, the observations appear to 
strongly favour a low density Universe, but the situation 
is more complicated if the galaxy luminosity distribution is 
strongly biased relative to the fluctuations in the density 
field. 

An extreme example of a biased model is the merging 
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model of Matarrese et al. (1997), in which visible galaxies are 
identified with very massive (~ 10 13 h^ 1 Mq) dark matter 
halos, with a linear bias relative to the mass fluctuations of 
b = 1.46 at z = and evolving as b = 0.41 + 1.05(l + z) 18 . As 
£(r) oc b 2 , this evolution of biasing would from z = to z — 1 
increase £(r) by a factor 7.76, producing clustering evolution 
of approximately e ga i axies ~ e mass — 3. Not suprisingly, this 
model overpredicts the observed faint galaxy co(9) and could 
be rejected in any cosmology (Moscardini et al. 1998). 

Matarrese et al. (1997) also present a 'transient model' 
in which star-forming faint galaxies are identified with less 
massive (~ 10 11 /i"'M ) halos and are assumed to evolve 
into a weakly clustered (b — 0.67 hence ro — 2 h _1 Mpc), 
low surface brightness population at the present day (see 
also Efstathiou 1995). The bias again evolves rapidly, as 
b = 0.41 + 0.26(1 + z) 1 ' 85 , which at z ~ 1 gives a similar 
£(r) to a model with normally clustered (b ~ 1) galaxies 
and e ~ 0. For a z < 1.6 sample (approximately the depth 
of our data) the transient model oj(9) is indistinguishable 
from an unbiased (b = 1) model, so our uj(6) results alone 
may not be sufficient to exclude the transient model. Higher- 
order correlation functions (Section 9.3) may provide further 
constraints on galaxy biasing. 

Limits on clustering evolution from this data will only 
apply to z ~ 1.5, as few galaxies will be detected beyond 
this redshift (Figure 11). The clustering of Lyman break 
galaxies at z ~ 3 (Giavalisco et al. 1998) is even stronger 
than expected from our e = model, indicating a negative 
e and an increase in linear bias at these redshifts. However, 
at least for the /-limited HDF sample, this increase in bias 
appears to affect u)(6) only at z > 2.4 (Magliocchetti and 
Maddox 1998), in contrast with the Matarrese et al. (1997) 
merging model in which bias rapidly increases at all z > 0. 
As only ~ 14 per cent of galaxies at the HDF limit are at 
z > 2.4, the u>(9) for the full /-limited sample (Villumsen et 
al. 1997) remains very low. 

Villumsen et al. (1997) originally interpreted the low 
ui(9) amplitudes on the Hubble Deep Field as favouring 
ro m 4 h^ 1 Mpc locally with e = 0.8 clustering evolution, 
whereas our results are much more consistent with e — 0. 
Furthermore, faintward of R ~ 26, it is unlikely that N(z) 
will become much more extended, and the most important 
change in N(z) is that less luminous galaxies are seen in the 
same volume of space. Hence it seems more likely that any 
further decrease in uj(9) amplitude at R > 26 is the result 
of a decrease in the intrinsic strength of clustering (ro) at 
low luminosities rather than clustering evolution. This in- 
terpretation would be supported by the results of Connolly, 
Szalay and Brunner (1998), who for photometric-redshift- 
divided HDF galaxies at / < 27 found weak intrinsic clus- 
tering of ro = 2.37 hT 1 Mpc but approximately stable clus- 
tering (e = —O^l^gg) over the 0.4 < 2 < 1.6 range (where 
dwarfs will dominate the sample) . 

We find no evidence of a flattening of the slope of uj(9) 
on going faintward (Table 4). None of our calculated to(9) 
are significantly flatter than 5 — 0.8, at any magnitude limit, 
which suggests that the £ (r) of the spiral galaxies forming 
the greater part of our catalog is close to 7 = 1.8, and cer- 
tainly no flatter than the 7 = 1.72 of LMEP. Giavalisco et 
al. (1998) estimate a S = 0.98±£ig for the u{0) of z ~ 3 
galaxies, suggesting the relatively steep slope we observe is 
maintained to much higher redshifts. 



A constant slope for to(9) may not be surprising in view 
of the Colin et al. (1997) simulations which predict no signif- 
icant flattening of the mass distribution £(r) , even to z ~ 5, 
for any Q. They do predict some flattening (A(7) ~ —0.2 to 
z ~ 1) of the galaxy £(r) for a strongly biased model, but 
such models already seem disfavoured by the u>(8) scaling. 
Some flattening of to{9) between R = 21 and R = 23.5 might 
be expected from the decrease in the fraction of ellipticals 
(Figure 12b), but we estimate that this would only amount 
to A(5) ~ -0.02 to -0.04 (from the LMEP and Guzzo et al. 
7 respectively), and might be cancelled out by the increase 
in the fraction of lower luminosity latte-type galaxies, as 
these may also have a steeper £(r) than L* spirals (LMEP 
estimate 7 = 2.01 ± 0.1). 

Further investigation of the evolution of u>(6) will re- 
quire large-format CCD surveys of a similar area to a greater 
depth, with imaging in three or more passbands (preferably 
including a near-infra-red passband) to enable a division of 
all the galaxies by both Hubble type and photometrically 
estimated redshift. 



9.2 Merger Rate Evolution 

On the larger of our two fields, we find a ~ 4cr significance 
excess in uj(6) at 2 < 9 < 5 arcsec compared to the (9~ 0,8 
power-law fitted at 9 > 5 arcsec. This was interpreted in 
terms of the fraction of the total sample of galaxies within 
9 < 5 arcsec pairs in excess of the expectation from the clus- 
tering seen at larger separations, which appeared to increase 
on going faintwards. The evolution in the pair fraction will 
follow approximately the evolution of the merger rate, pa- 
rameterized here as R m = R m o(l + z) m . 

By comparing the fraction of galaxies in close pairs at 
V < 22.5 with that at very low redshifts, Carlberg et al. 
(1994) claimed rapid evolution of m = 3.4 ± 1.0, but this 
was based on a very small sample. With a dataset ~ 40 
times larger, Infante et al. (1996) estimated a merger rate 
evolution of m = 2.2 ± 0.5 at R < 21.5, and we similarly 
estimate m = 2.0lto'gg f° r our field ' e ' data to R — 23.5. 

The Infante et al. (1995) results and ours are consistent 
with the m ~ 2 expected for a low density (ft ~ 0.2) Uni- 
verse, but inconsistent with the m ~ 4 predicted for f2 = 1 
(Carlberg et al. 1994). Neuschaefer et al. (1997) derived an 
even lower rate of m = 1.2 ± 0.4 from the close pair frac- 
tion on deeper (/ < 25) HST data, which may hint at a 
reduction in m with increasing survey depth, and would be 
consistent with the Carlberg et al. (1994) model for a low £1 
Universe in which merging leads to a moderate reduction in 
mean galaxy mass at high redshifts. 



9.3 Hierarchical Moments and Non-linear Bias 

Using a simple counts-in-cells method, we investigated the 
area-averaged three and four point correlation functions of 
18.5 < R < 23.5 galaxies, in the form of the hierarchical 
moments S3 (9) and 34(8). Both moments were detected with 
high statistical significance, and, over a range of cell sizes 
0.01 < 9 < 0.1 deg, showed relatively steep slopes best- fitted 
with power-laws 9~ s with 5 = 0.413 ± 0.083 for s 3 (9) and 
S = 0.614 ± 0.134 for sa(9). These slopes appear consistent 
with those of the EDSGCS sj(8) at equivalent scales (0.15 < 
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r < 1.5/i _1 Mpc), but strongly inconsistent with the near 
zero or positive 8 of the APM sj(9). 

Szapudi and Gaztanaga (1998) attribute the difference 
between the APM and EDSGCS sj(8) at small scales to 
the more efficient deblending of galaxy images in rich clus- 
ter fields in the EDSGCS analysis, but also claim that the 
EDSGCS contained some spurious detections, and so con- 
clude that 'the true galaxy distribution is likely to lie some- 
where between the two surveys'. Our results, based on CCD 
data giving much more reliable deblending than digitized 
photographic plates, suggest the EDSGCS results are the 
more accurate. 

Furthermore, a recent direct measurement of the three- 
point correlation function of 0.03 < z < 0.15 galaxies (Jing 
and Borner 1998) in the Las Campanas Redshift Survey 
(also selected using CCD data) found the ratio of the pro- 
jected three-point correlation function to £(r) to scale with 
separation as r~ 0,3 at 0.2 < r < 3/i _1 . These results were 
consistent with the slope and normalization of the EDSGCS 
83(8), and with our best-fit slope for 83(6), but are inconsis- 
tent with the flat APM s 3 (9). 

The Sj(r) from iV-body simulations (Baugh, 
Gaztanaga and Efstathiou 1995) do show some rise towards 
smaller cell sizes, which is steeper for successively higher or- 
der moments, as observed here and by Szapudi et al. (1996). 
The simulations give no suggestion that the Sj(r) would 
flatten between z mean — 0.69 and z mean — 0.14 epochs, but 
predict a slight steepening with time, and so may strengthen 
the case that our results favour the EDSGCS sj{9) over the 
APM estimate. The predicted Sj(r) are steeper for lower 
density cosmologies (Gaztanaga and Baugh 1995); our re- 
sults are more consistent with the SI = 0.2 than the SI = 1 
models, as is the three-point correlation function of Jing and 
Borner (1998). 

At this point many of our results (e ~ galaxy clus- 
tering, the moderate evolution of the merger rate, the steep 
•83(0) and S4,(9)) seem to favour a low SI cosmology, but 
this interpretation assumes that galaxies at least approxi- 
mately trace the underlying mass distribution. Strongly bi- 
ased models (e.g. Matarrese et al. 1996; Moscardini et al. 
1998) might be adjusted to fit the ui(9) scaling for a wide 
range of cosmological models, but, as we discuss below, 
higher-order correlation functions may have the potential 
to break the apparent degeneracy between the biasing and 
cosmological model. 

For a distribution of mass p ma ss(r), with fluctuations 
about the mean <5 mass (r) = p mass (r) — (p mass ), and a dis- 
tribution of the galaxy luminosity in a particular passband 
Pgal(r), with fluctuations S ga i(r) = p ga l{r) - (p ga i), the bi- 
asing of galaxies relative to mass can be expressed as the 
series (e.g Fry and Gaztanaga 1993) 

S ga i(r) = J2^ 5 rnass(r) (27) 

fc=0 

The galaxy £(r) can be interpreted as the £(r) of the un- 
derlying mass distribution multiplied by b\, where b\ is the 
linear bias. The hierarchical moments of the galaxy distri- 
bution are related to those of the mass distribution through 
expressions involving the higher-order (non-linear) bias co- 
efficients, 



s 3(gal) = (S3(mass) + 3 — )/6l (28) 

b b b^ 

S 4(gal) = (S 4 f ma3s \ + 12 — S 3 ( mass ) + 4 — + 12-2 )/^l (29) 
01 Ol o 1 

Gaztanaga and Frieman (1994) found the APM survey hier- 
archical moments to be consistent with perturbation theory 
predictions for a simple model in which the blue-band lu- 
minosity from galaxies linearly traces the mass distribution 
- i.e. 61 ~ 1 and bk — for all k > 1 - with the best fit 
for a low density (SI = 0.2 A = 0.8) CDM model. Improved 
modelling with TV-body simulations (Baugh, Gaztanaga and 
Efstathiou 1995; Gaztanaga and Baugh 1995) gave steeper 
Sj(r) at r < 7/t -1 Mpc for a given SI, but appears to give 
the same interpretation when compared with the EDSGCS 
results (Szapudi and Gaztanaga 1998). 

However, even if an optically-selected sample as a 
whole is approximately unbiased, the different types of 
galaxy may vary greatly in biasing properties. For exam- 
ple, 60pm-selected IRAS galaxies, in addition to having a 
lower £(r) amplitude than optically selected galaxies (im- 
plying &i(/ijAS) — 0.7), are a factor of 2 lower in 83(8) and a 
factor 3 lower in S4(9) (Fry and Gaztanaga 1993). This indi- 
cates that IRAS galaxies have a non-linear bias with b 2 < 0, 
so that they tend to avoid the richer clusters in the galaxy 
distribution as a whole. 

In simulations (Baugh, Gaztanaga and Efstathiou 1995; 
Gaztanaga and Baugh 1995), the Sj(r) of the underlying 
mass distribution slowly increase and steepen with time at 
the r > 7/i~ 1 Mpc distances relevant to our survey. At a 
fixed proper separation r ~ Ih^ 1 Mpc, the models appear 
to show a decrease of ~ 20-25 per cent for Ss(r) and ~ 30- 
40 per cent for Si(r) on going from the mean epoch of the 
EDSGCS/APM galaxies to that of the 18.5 < R < 23.5 
data. The decrease we observe is consistent with this model, 
implying that the mean linear and non-linear biasing prop- 
erties of the luminosity from galaxies (in the rest-frame blue- 
band) galaxy remain approximately constant to z ~ 1. 

We interpret this as evidence that (i) most rich galaxy 
clusters formed prior to our mean redshift 2 ~ 0.7, as other- 
wise the galaxy distribution at these redshifts would be more 
like local IRAS galaxies, resulting in a further reduction of 
as much as a factor of ~ 2-3 in the sj(9) and (ii) galax- 
ies with different bias properties (cluster and field galaxies, 
ellipticals and spirals, dwarf and giant galaxies) undergo a 
similar evolutionary brightening in the rest-frame blue-band 
out to z ~ 1, as in PLE models. 

Our results would argue against 'transient' or 'disap- 
pearing dwarf models (e.g. Babul and Rees 1992; Efstathiou 
1995), in which a large proportion of the faint blue galax- 
ies are weakly clustered low mass systems undergoing brief 
starbursts. Although the 'transient model' of Matarrese et 
al. (1996) would be consistent with the u(ff) at R = 23.5, 
the transient dwarfs are associated with low mass haloes col- 
lapsing at relatively late epochs, which would have a very 
negative 62 (Mo and White 1996). Domination of the deep 
sample by such objects would therefore change the nonlin- 
ear biasing properties relative to the APM/EDSGCS galax- 
ies and give 83(9) and S4,(6) lower than are observed. PLE 
type evolution is also favoured by combined HST and spec- 
troscopic data which provide direct evidence that field and 
cluster galaxies (Schade et al. 1996a, 1996b) and early and 
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late types (Roche et al. 1998; Glazebrook et al. 1998) all 
undergo similar blue-band evolution of ~ 1 mag to z ~ 1, 
so that, at least in visible-light wavelengths, the evolution 
required to fit the number counts is largely shared amongst 
normal galaxies rather than being concentrated into a new 
type of object. 

Furthermore, if galaxies, on average, linearly trace the 
mass distribution at both low redshifts and the depths of 
our survey, the faint galaxy uj(0) scaling indicates that clus- 
tering is stable (e ~ 0) for the mass distribution and not 
just for the galaxy luminosity as observed in a particular 
passband. Hence the hierarchical moments strengthen our 
earlier conclusions that galaxy clustering is genuinely stable 
to z ~ 1 and that we live in a low f2 Universe. 

As larger areas are surveyed at CCD depths, the hier- 
archical moments are likely to become an important com- 
plement to the study of the faint galaxy uj{9). We hope to 
further these investigations, As for lo{6), by surveying a sim- 
ilar area to a greater depth and by observing in three or 
more passbands, including one in the near-infra-red, so that 
colour-divided samples will better constrain the evolution 
of sj{6) and show the relation between non-linear bias and 
galaxy colour. 
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